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We study the two-dimensional generalized Weingarten model reduced to a point, which 
interpolates between the reduced Weingarten model and the large-A^ gauge theory. We 
calculate the expectation value of the Wilson loop using a Monte Carlo method and determine 
the string tension and string susceptibility. The numerical result suggests that the string 
susceptibility approaches —2 in a certain parametric region, which implies that the branched- 
polymer configurations are suppressed. 



§1. Introduction 

In this paper, we study a model defined by the action 

s = -? N E Tr (44v^) + aN E ^ ( A l A » + ' ( M ) 

where (/x = 1,2, ••■ , d) are complex N x N matrices. (In t he a ctual numerical calculation 
describe in § [31 we consider only the case d = 2.) This mode l***)! interpolates 1 ) between the 
reduced Weingarten model 3 ) and the large- N reduced U(N) gauge theory. 4 ) We regard the model 
(ll-ip as describing an ensemble of random surfaces and study whether it can describe the Nambu- 
Goto string. 

The original Weingarten model 5 ) was proposed as a nonperturbative description of the Nambu- 
Goto string. This model is defined as follows. Consider a d-dimensional square lattice Z d and 
introduce a complex N x N matrix A x ^ for each link connecting the sites x and x + fx in such a 
way that A x+ f l ^ )1 = A X:fl . Then the action of the Weingarten model is given by 

d 

s w = -npJ2Y,^( a ^ a *+^ a U>,A,v) + iV EE Tr (4/V^) • (i-4) 

X X fj,= l 



*' E-mail: hanada@gauge.scphys.kyoto-u.ac.jp 
**' E-mail: kubo@gauge.scphys.kyoto-u.ac.jp 

***) This model is the zero-dimensional reduced version of the "interpolating model" proposed in Ref. 1). For 
a > 1, this model is equivalent to that studied in Ref. 2), 

d d 2 

S = -7iV^Tr(xtAlX M X„) + K iV ]T Tr - l) , (1-2) 

(J-J^V jit— 1 

with the relation 

A » = f-l x »> ft = a { 1 -^) 2 ' ^(T^lf (1 ' 3) 
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The partition function is given by 

Z\y = J drriN exp ( — Sw) , (1-5) 



where the measure dmN is defined by 



dm N = II II (7 d [ Re (4r, M )ij] d pmfAr^y] j . (1-6) 

Let Cj represent closed contours on the lattice. Then, multiplying A M along d and taking the 
trace, we obtain Wilson loops w(Ci). The correlator of w(Ci), defined by 

W(d, ■■■ ,C n ) = ^- J dm N exp (S W ) ^w(C{) ■ ■ ■ ^w(C n ), (1-7) 

is evaluated as 

W(d,--- ,C n ) = N 2 ~ 2n exp(-a( S )log/3- 1 -Ms)logiV 2 ), (1-8) 

s6S({C0) 

where S({Ci}) is the set of surfaces on the lattice whose boundary is C\ U • • • U C n , a(s) is the 
area of the surface s, and h(s) is the number of handles of s. If we regard log/3 -1 and as the 
string tension and the string coupling, respectively, then Eq. (j 1 • 8[) can be interpreted as the sum 
of random surfaces weighted by the Nambu-Goto action. 

Next, let us consider the reduced Weingarten model, 3 ) whose action is given by 

d 

Srw = -N0 ]T TV (^A„44) + N £ Tr . (1-9) 

This action is invariant under the U(l) d transformation 

- e*M M . (MO) 

If this symmetry is not broken spontaneously in the large- iV limit, the correlators of the Wilson 
loops in this model are identical to those of the original Weingarten model, (|l-4p . Because the 
reduced Weingarten model has only d matrices, numerical calculations are more tractable. This 
model was studied numerically 6 ) in the cases d = 2,3, and it was shown that the Weingarten 
model does not describe smooth surfaces but, rather, branched polymers. 7 ) 

One possibility to overcome this difficulty is to consider the modified action (jl-ip . This action 
is motivated by the following observations. First, this model interpolates 1 ) between the original 
Weingarten model (a = )*'\ and the reduced U(N) gauge theory 4 ) (a = oo). Because both of 
those models are thought to be related to string theory, it is natural to consider the intermediate 
region. In this region, this model allows a lattice string interpretation similar to that of the original 
Weingarten model, because the relation (jl-8p holds also in this model, as long as the surface s does 



*- ) We need the redefinitions 



^(original) _ ^^(modified) ^(original) _ 1 ^(modified) ^.-^ 
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not intersect itself. Second, the action (jl-ip becomes a set of d copies of a complex one-matrix 
model with a double-well potential. In the case of the Hermitian matrix model, we can describe 
a type OB string by flipping the sign of the double-well potential. 8 ) Therefore we conjecture that 
also in the case of the Weingarten model, a worldsheet supersymmetry can be introduced by 
modifying the potential, and this may prevent the worldsheet from becoming branched polymer. 

This model has been solved analytically only in the special cases P = 1 ) and a = oo for 
d = 2. 9 )' 10 ) In Ref. 2), the parametric region a > 1 is studied. For sufficiently large a, there are 
d phase transitions that correspond to the partial breakdowns of U(l) d symmetry. For d > 3, 
these phase transitions smoothly approach the known phase transitions of the large- N reduced 
U(N) gauge theory 11 )' 12 ) as a — > oo. For sufficiently small a (a < 2 in the case d = 2), these two 
transitions seem to occur simultaneously. In this paper, we study the parametric region a < 1 in 
detail in the case d = 2. For a technical reason, we study the maximally twisted model. 

The organization of this paper is as follows. In § [2j we present theoretical preliminaries. In 
§EJ we present the numerical result. We study the phase diagram in § 13.11 and then determine the 
string tension and string susceptibility in § 13.21 Section [4] is devoted to conclusions and discussion 
of future directions. In Appendix [A] we comment on numerical results concerning the generalized 
Weingarten model without a twist. 

§2. Theoretical preliminaries 

2.1. Parametric region of the generalized Weingarten model 

We begin by considering the parametric region in which the model defined by the action (|1-1D 
is well-defined. As shown in REf., 2 ) the action is bounded from below if and only if 

£ ih- (24 » 

This can be seen as follows. In order to determine whether the action is bounded, it is sufficient 
to consider the quartic term, 

d d 

SUartic = ~PN ]T Tr (AlAlA^Au) + aN ^ Tr (4^4^) . (2-2) 
Then, using the inequality 

2Re Tr (AB^ < Tr (aaA + Tr (^-B f ) , (2-3) 

we have 

2Re Tr (A^A^Ai) < Tr (a^A^AI) + 

< Tr (4^4^) + (AlA v AlA v ) . (2-4) 

Summing over the spacetime subscripts, we obtain 

d d 2 

^ (44 < ( d ~ !) E Tr (44) • (2-5) 
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Combining this relation with (|2-2p . we obtain the bound (|2-ip . This bound is indeed realized in 
the case of the unit matrix, for example. 

For [3a~ l > -py, although the action is not bounded from below, the model still can be 
well-defined for large-iV. 3 ) In the original Weingarten model, the free energy per unit volume, 
F(N,{3), is given by 

1 oo OO / A \ 

g=0 A=0 

where n g (A) represents the number of closed surfaces with genus g and area A. By taking the 
planar limit N — > oo, with (3 kept fixed, we obtain 

1 oo 

—F(N,P) = Y,MA)f3 A . (2-7) 

In Refs. 13) and 14), it is argued that n${A) behaves as 

n (A) ~ A h ~ x k A (2-8) 

for sufficiently large A, where b and k are universal and regularization-dependent constants, re- 
spectively. Then, we have 

-. oo 

^F(iV,/3)~ J^ 6 " 1 ^, ( 2 - 9 ) 



A=0 



and this quantity is finite for (3 < (3 C = p In the same way, it can be shown that the expectation 
value of the Wilson loop is also finite. Therefore, the original Weingarten model is well-defined 
for [3 < 0c in the large- iV limit. 3 ) 

In numerical simulations, obviously we cannot make iV infinite. However, for large enough N, 
there is a metastable state corresponding to the planar limit. The "lifetime" of this metastable 
state becomes longer as N increases. 6 ) 

For generic a, the measure in the strong coupling expansion changes from the Gaussian one 
in the original Weingarten model as 



exp ^-Njy^ A l A ^j dm N — > ex P ^-aNj2^ (^A " 1 + ^) j 



dm N . (2-10) 



Because we cannot evaluate the strong coupling expansion exactly, we simply assume that the 
effect of the change of the measure can be absorbed into changes of b and (3 C : 

b, (3 C — ► b(a), f3 c (a). (2-11) 

If this assumption is correct, then the planar limit exists for [3 < (3 c (a) also in this case. As we 
see in § [3l the numerical data seem to be consistent with this assumption. 

2.2. String tension and string susceptibility 

First, let us consider the case of the original Weingarten model. Then, the number of planar 
random surfaces on a lattice with boundary C and area A is given 

by 13),14) 

n (A; C) ~ A b k A = A b f3~ A (2-12) 
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for large A. Note that no (A; C) is A times larger than no(A), because there is a degree of freedom 
corresponding to the choice of the location of a puncture. Then, the expectation value of the 
Wilson loop is evaluated as 

W(C) = J2MA;C)(3 A 

~ const, x Y A \P/Pc) A 
A>A 

~ W C {C) - const, x |/3 - () c \~ b ~ l + • • • , (2-13) 

where W C {C) = W(C)\p = p c and Aq represents the minimum area surrounded by C. Because the 
relation (|2-12p can hold only for sufficiently large A, only the leading order singularity in the limit 
(3 — ► p c is reliable in the expression ()2- 13|) . 

We can readily determine b from (|2T3p . In Refs. 15) and 6), b is determined to be —1.5. 
However, as we see later in this section, branched polymers dominate the path integral if b > — 2. 7 ) 
Furthermore, the string tension is finite even at j3 = /3 C , 15 ^ although, in order to take the continuum 
limit, the string tension must approach zero. For these reasons, the original Weingarten model 
does not allow a meaningful continuum limit. 

Next let us consider the parametric region a > 0. In this case, if we assume the number of 
random surfaces is changed effectively as Eq. (|2-lip . then we can expect the following behavior at 
(3 ~ (3 c {a): 

W{C)~ Y, ^ a Xf3/f3 c (a)) A - 

A>A 

~ W C (C) - const. x\/3- [3 c (a)\- h(a) - 1 + ■■■ . (2-14) 

In §[3] we determine the value of the "string susceptibility", b(a), on the basis of this relation. 

In the latter part of this section, following Ref. 7) 
we show that smooth surfaces (resp., branched poly- 
mers) dominate the path integral if the string sus- 
ceptibility b is smaller (resp., larger) than —2. The 
number of planar surfaces with sufficiently large area 
A is given by Eq. f|2-8j) . Now, let us assume that Fig. l. A surface with a pinch, 

smooth surfaces are dominant for a given b. Then, 
the number of surfaces with a pinch (Fig. [T|) is 

A A 

Y MA'; C)n {A - A'; C) ~ k A Y A ' b ( A ~ A 't ~ k A A 2b+1 , (2-15) 
A'=0 A'=0 

where C represents the punctures to be connected. In order for the smooth surfaces to dominate 
the pathintegral, the number of smooth surfaces no (A) must be larger than the number of surfaces 
with a pinch. Therefore, we have the following relation: 

n (A) ~ k A A b ~ l > k A A 2b+l ^b<-2. (2-16) 

By contrast, if b > —2, surfaces with more pinches contribute more, and hence branched polymers 
dominate the path integral. 
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§3. Numerical results for the two-dimensional generalized Weingarten model with 

maximal twist 

In this section, we present the numerical result for the two-dimensional generalized Weingarten 
model. In order to determine the string susceptibility using the ansatz (|2-14j) . we need to study 
the metastable region. Therefore, in order to suppress the tunneling effect, 6 ) we studied the 
maximally twisted reduced model. 16 ) This is obtained by making the replacements 

0^-0, W[m,n] -» (-) mn W[m,n], (3-1) 

where W[m, n\ represents the expectation value of an m x n rectangular Wilson loop. 
3.1. Phase diagram 

A schematic picture of the phase diagram for the maximally twisted reduced model is displayed 
in Fig. [2 The line = C represents the boundary of the metastable region; i.e., for < j3 c , 
the "lifetime" of the metastable state becomes longer as N becomes larger. As the value of a 
increases, this line approaches the boundary of the stable region, a = 0. For a ~ 1.2, a phase 
transition takes place at = /^breakdown- The values of /^breakdown seem to diverge as a — ► oo. This 
is consistent with the analytic result at a = oo (i.e., the unitary limit), according to which no 
first-order or second-order phase transition exists. 

A remark is in order here. As shown in Appendix [AJ ^breakdown seems to coincide with 0\ 
in untwisted model. This represents the breakdown of the U(l) 2 symmetry. Furthermore, the 
expectation value of the Wilson loop seems to coincide with that in the untwisted model not only 
for < /3brcakdown = Pi but also for > /^breakdown = 01- Therefore, it is plausible that the 
phase transition at = /3breakdown corresponds to the breakdown of the U(l) 2 symmetry. If this 
is indeed the case, then the large- N reduction 3 )' 4 ) cannot be applied for > 0breakdown and hence 
the lattice-string interpretation would not be possible for this model in this parametric region. 




first order? 

Fig. 2. Phase diagram of the two-dimensional generalized Weingarten model with a twist for a < 1. 
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Fig. 3. W[l, 1] for a = 3.73, /3 = 0.4. The horizontal 
line represents the estimated value in the large- 
N limit, obtained by fitting the form W[l, 1] = 
ci + C2N C3 , where ci,C2 and C3 are constants. We 
find good convergence at N = fOO. 
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Fig. 4. W[5, 5] for a = 3.73, /3 = 0.4. The horizontal 
line represents the estimated value in the large- 
N limit, obtained by fitting the form W[5, 5] = 
ci + C2N° 3 , where ci, C2 and C3 are constants. We 
cannot find good convergence at this level. 



3.2. Wilson loop 

Let W[m, n] be the expectation value of the m x n rectangular Wilson loop. We calculated 
the values of W[m, n] with m,n = 1, • • • ,5. For each a and (3, we fitted the data numerically by 
using the ansatz 

W[m, n] = c- p 2 ^ m+n ) ■ e - T - mn , (3-2) 

where c,p and T are constants. There is a subtlety here: The numerical result suggests that the 
effect of the perimeter p depends not only on a but also on j3. Therefore, we expect the behavior 
described by Eq. (|2-14j) only when the loop is large enough. In numerical simulations, because we 
can only study small loops, we should take the perimeter effect into account. Then, we expect 

p~ L W[m,n] = W'[m,n] ~ W' c [m,n] - const, x \(5 - p c {a)\~ b{a) ~ 1 + ■■■ (3-3) 

instead of Eq. (|2T4p . As we see below, this ansatz appears to be consistent with the numerical 
data. 

For N = 100, only the loops satisfying mn < 10 seem to converge (see Figs. [3] and [4]). For 
this reason, we use only such values to extract the string tension and the string susceptibility. 

3.2.1. a = 0: Original Weingarten model 

For a = 0, it is known that 6 )' 15 ) 

b =-1.5. (theoretical) (3-4) 

Then, taking the perimeter effect into account, we find that the numerical data are consistent 
with this value (Fig. [5]). From the value of W[l, 1], we obtained 



b = -1.50 ± 0.02. (numerical) 



(3-5) 
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Fig. 5. log |Wo[w., n] — W[m, n]\ as a function of log |/3 — (3 C \ for q = (original Weingarten model), with N = 100. 
Here, the perimeter effect is taken into account. The slope of the line is 0.5. [Left] (m,n) = (1, 1). [Right] 
(m,n) = (2,1). 
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log W[m, n] as a function of the A — mn for a = 1,N — 100. [Left] Without the perimeter effect taken 
consideration. [Right] With the perimeter effect taken into consideration. The slope corresponds to — T. 
error bars become larger due to ambiguity in p. 



3.2.2. 0<a<1.2 
String tension 

Taking the perimeter effect into account, we can determine the string tension T (see Fig. [6]). For 
fixed a, the string tension T decreases as /? — > f3 c (ot) (Fig. [7]). Although the value T (a, f3 c (a)) 
decreases as a becomes large, it remains positive (see Fig. [5]). 
String susceptibility 

We next determine the "string susceptibility" b(a) using the numerical data and the expression 
(|3 3f) . As can be seen from Fig. [H the ansatz (]3-3|) agrees better with the data than the ansatz 
(|2-14p . Figure [TUl displays a plot of log \W C [1, 1] - W[l, 1]| as a function of log |/3 - (3 C \ for a = 1. 
We can read off the value —6 — 1 from the slope of the fitting line. The value of b{a) is plotted 
in Fig. [12j For a ~ 1.2, 6(a) seems to approach —2 (Fig. [T2]) . The uncertainties here are rather 
large, and they come mainly from an ambiguity in the perimeter effect p. One of the reasons that 
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Fig. 7. The string tension T for a = 1.2, N = 100. 
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Fig. 8. The string tension T for f3 = c (a),N = 100 
as a function of a. 
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Fig. 9. W c [l, 1] — W[l, 1] is as a function of f3 — P for a = 1, N = 100. The curves represent the fitted form 
W = W c — const • |/3 — f3 c \^ h ^ 1 . [Left] With the perimeter effect taken into consideration. [Right] Without the 
perimeter effect taken into consideration. 



they are so large is that the deviation of the value of the Wilson loop (at N = 100) from that in 
the large- N limit is not uniform but depends on the size of the loop and the parameters a and (5. 
If we can use larger matrix size N, we would be able to make the uncertainties smaller. 

3.2.3. 1.2<a<oo 

In this parametric region, a phase transition takes place at (3 = /^breakdown < flc- As remarked 
in § 13-H the U(l) 2 symmetry seems to be broken for > /3breakdown- Therefore, it is plausible 
that the reduced model deviates from that on a lattice, and thus the lattice-string interpretation 
is not possible!*^ 



Similar phenomena in the unitary gauge theory are studied in Ref. 17). 
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Fig. 10. log|W c [l,l] - W{1,1)\ as a function of 
log (5 C | for a = 1, N = 100. Here, the perimeter 
effect is taken into consideration. The line repre- 
sents the fitted form, whose slope is — b — 1. 
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Fig. 11. W[l, 1] as a function of (3 for a = 1.2, JV = 
100. Here, the perimeter effect is taken into ac- 
count. 
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Fig. 12. String susceptibility b as a function of a. The result, 6 
result found in Refs. 15), and 6). 



-1.50 ± 0.02 at a — 0, is consistent with the 



§4. Conclusions and discussions 



In this paper, we have reported the results of a a numerical study of the two-dimensional 
generalized Weingarten model (|l-ip . If we assume the relations (j2- lip and (j3 3h . then the numerical 
data show that the string susceptibility approaches —2 as we increase the value of the parameter 
a to ~ 1.2. This result suggests that branched-polymer configurations are suppressed in this 
parametric region. We also found that the string tension decreases. However, we encountered a 
phase transition before the string tension becomes zero. Therefore, we cannot take the continuum 
limit in the reduced model studied in this paper. 

There are several pontential directions for future studies. First, it is necessary to understand 
the present model in terms of the random surface and explain why the assumption (|2-lip seems 
to be consistent with the numerical data. Doing so, we would be able to clarify whether or not 
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our numerical result actually indicates that the branched polymer configurations are suppressed. 
In addition, such an understanding would be helpful for finding better models. Second, it would 
be interesting to study the model defined on a lattice. ^ In the case of the large- N reduced d- 
dimensional U(N) gauge theory (d > 3), the breakdown of U(l) d is an artifact of the reduced 
model; if the model is defined on a lattice of size L d , then the U(l) d symmetry is not broken to 
the weaker coupling as L increases. 11 ) If similar phenomena exist in the present case, then using 
the model with L > 1, we could study the larger parametric region, and we may be able to find a 
point at which we can take the continuum limit. Third, we can also consider higher-dimensional 
models. In this case, the lattice-string interpretation may be valid in the parametric region in 
which some of the U(l)s remain unbroken. This point is worth studying. We hope to report 
analysis of these models in future publications. 
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Appendix A 

Comparison with the Two- Dimensional Generalized Weingarten Model without a Twist 



A.l. Two-dimensional generalized Weingarten model without a twist 




{1}^ [7(1)2 

Fig. 14. Phase diagram of the two-dimensional gener- 
Fig. 13. Phase diagram of the two-dimensional gener- alized Weingarten model without a twist for a > 1. 
alized Weingarten model without a twist for a < 1. (This figure is based on Fig. 10 of Ref. 2).) 
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without a twist y 
with a twist t- 



without a twist , — | — | 
s s with a twist \ — x — s 



Fig. 15. Plot of the expectation value of 1 x 1 Wilson Fig. 16. Plot of the expectation value of the action 
loop for a' 1 = 0.3, N = 50. for or 1 = 0.3, N = 50. 



In this section, we consider the original, untwisted generalized Weingarten model in two 
dimensions. 

For large, fixed a there are two curves of first-order phase transitions. We call them (3\ and Pi 
in ascending order. They correspond to the breakdown of the U(l) 2 symmetry. If we increase (5 
with a fixed, U(l) 2 is broken to U(l) at f3\, and then it is broken completely at Pi- The values of 
Pi and Pi seem to diverge as a -1 — ► 0. This is consistent with the analytic result, in which U(l) 2 
is not broken. At small a, (5\ and Pi seem to become equal. This transition persists to a ~ 1.2 
where it merges with the boundary of the metastable parametric region. Near the boundary of 
the well-defined region f3 < a, the U(l) 2 symmetry is restored. In Fig. [T3l this line of restoration 
of U(l) 2 is drawn by a dotted line near a ~ 1.2, because the restoration cannot be seen clearly in 
this region. 

In the parametric region where the U(l) 2 symmetry is not broken, our reduced model is 
equivalent to the model defined on a lattice 1 ) through the large- N reduction. If this symmetry is 
broken, then the reduced model deviates from that on a lattice, and the lattice-string interpretation 
is not valid Rl. 

A. 2. Comparison of untwisted and twisted models 

In this subsection, in order to avoid confusion, we indicate the physical quantities in the 
twisted model with the subscript T. For example, Wilson loops are denoted as 

W[m,n](a,(3), (without twist) (AT) 
W T [m,n](a, /3) = (-) mn W[m, n](a, -/?). (with maximal twist) (A-2) 

Here we have included the phase (— ) mn so that the two prescriptions give the same value in the 
strong coupling region: 

W[m, n](a, P) = Wr[m, n](a, P). (strong coupling) (A-3) 

*' For d > 3, if U(l) d remains unbroken, then this model is equivalent to a model defined on a d'-dimensional 
lattice coupled to d — d! adjoint scalars. For this reason, such a parametric region is of interest in this case. Note 
that the U(l)s do not necessarily break one-by-one in the case of the twisted model. 
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Fig. 17. Plot of the specific heat for a' 1 = 0.3, N 
50. 



Fig. 18. Plot of the expectation value of £ |TrA M | 2 
for a' 1 = 0.3, iV = 50. 



In Fig. [151 w e plot the expectation values of 1 x 1 Wilson loops for a -1 = 0.3. We see that W[l, 1] 
and Wt[1, 1] indeed take the same value. We can also see that the expectation values of the action 

are nearly equal (Fig. ITS]) . 

It is interesting that the two prescriptions seem to give the same expectation values for the 
Wilson loops not only in the strong coupling region but also in the weak coupling region, which 
is separated from the strong coupling region by a phase transition; indeed, as can be seen from 
Fig. [T71 phase transitions take place naer (3\(a~ l = 0.3) = 1 both in the twisted and untwisted 
models. In the case of the untwisted model, this transition corresponds to the breakdown of 
U(l) 2 (see Fig. US). In the case of the twisted model, because the expectation value of (Tr^l 
remains nearly equal to zero, this transition does not correspond to the complete breakdown of 
one of the U(l)s; it may represent a breakdown to Z n for some integer nF^l 
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A similar conjecture is made in Ref. 17) in the case of the four- dimensional twisted Eguchi-Kawai model. 



